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Abstract 

We investigate two approaches to increase the efficiency of phenotypic 
prediction from genome-wide markers, which is a key step for genomic se- 
lection (GS) in plant and animal breeding. The first approach is feature 
selection based on Markov blankets, which provide a theoretically-sound 
framework for identifying non-informative markers. Fitting GS models us- 
ing only the informative markers results in simpler models, which may allow 
cost savings from reduced genotyping. We show that this is accompanied by 
no loss, and possibly a small gain, in predictive power for four GS mod- 
els: partial least squares (PLS), ridge regression, LASSO and elastic net. 
The second approach is the choice of kinship coefficients for genomic best 
linear unbiased prediction (GBLUP). We compare kinships based on differ- 
ent combinations of centring and scaling of marker genotypes, and a newly 
proposed kinship measure that adjusts for linkage disequilibrium (LD). 

We illustrate the use of both approaches and examine their performances 
using three real-world data sets from plant and animal genetics. We find 
that elastic net with feature selection and GBLUP using LD-adjusted kin- 
ships performed similarly well, and were the best-performing methods in 
our study. 

Keywords: genome-wide prediction, genomic selection, feature selec- 
tion, Markov blanket, linkage disequilibrium, kinship. 
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1 Introduction 



The ever-increasing amount of genetic information available in plant and animal 
breeding is reflected in the development of sophisticated models for the predic- 
tion of quantitative traits from genome-wide markers, also known as genomic 
selection (GS). The markers are typically dense single-nucleotide polymorphisms 
(SNP). Approaches to the problem have moved from models w ith simple spec- 
ifications, such as ridge regression (IHoerl and Kennardl Il970l) and the LASSO 



(ITibshirani, 119961), to models based on highly- structured hierarchical distribu- 



Guan and Stephens! (1201 11) . 



tions, as in 

This complexity is motivated by the need to correctly model the genetic ar- 
chitecture of the trait under investigation, while producing models that are easy 
to estimate even for large SNP profiles. We focus on two key aspects of these 
models: the inclusion of a preliminary step that removes SNPs that appear to 
be redundant, and the choice of kinship matrices to model the relatedness of the 
genotyped individuals. 

The former is equivalent to feature selection (IKoller and Sahamil . 1 19961) . and 
can be achi eved by shrinking SNP effects towards zero, either through the use of 
constraints (IZou and HastieL 120051) or thr ough appropriate prior distributions in a 
Bayesian setting (IMeuwissen et al.l . l200lh . We examine the effectiveness in GS of 
Markov blankets (|Pearll Il988|) . which have been extensively studied in graphical 
modelling. They provide a principled solution to feature selection problems, and 
can be implemented as a data pre-processing step prior to fitting the GS model. 
We implement Markov blanket feature selection within four GS models applied 
to three real-world data sets covering barley, rice and mouse genetics. 

Kinship matrices were traditionally derived from pedigrees using a single def- 
inition, but with kinship s now being calculated f rom SNP data many different 
definitions are available (|Astle and BaldingL |2009|) . We investigate four kinship 
matrices within genetic best linear unbia s ed pred iction (GBLUP). These include 
a novel matrix introduced by ISpeed et al.l (|2012a|) which adjusts for the bias intro- 
duced by differences in local linkage disequilibrium (LD), and has been shown to 
increase the precision of heritability estimates. 



2 Background 



2.1 Markov Blankets and Feature Selection 

The Markov Blanket of a variable of interest T, denoted as B(T), is the minimal 
set of variables conditioned on which al l other variab les in the model are prob- 
abilistically independent of the target T (IPearll . 1 19881) . The Markov blanket of a 
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phenotype y in a GS model is the minimal set B(y) C X such that 



P(y|X)=P(y|5(y)), 



(1) 



that is, the subset of SNPs B(y) that makes all other SNPs redundant as far as the 
trait y is concerned. Given this property, knowledge of only the SNPs in B(y) is 
enough to determine the probability distribution of y. Other SNPs become super- 
fluous, either because they are not associated with the trait or because their effect 
is mediated by the SNPs in B(y). If B(y) were known, any GS model could be 
fitted using B(y) instead of the full SNP profile X with no loss of information, but 
in practice the need to estimate B(y) means that some informati on loss is possible . 



This t wo-stage approach contrasts with mo dels such as BayesB (fMeuwissen et al. 



200 lh and the LASSO (lTibshiranii 119961) . which select significant SNP effects 
concurrently with model fitting and in a model- specific way. 

Markov blankets can be estimated efficiently from data through the use of con- 
ditional ind ependence tests, such as parametri c and non-parametric te sts for partial 
correl ation (fLegendrel BOOOl iHotellingl 1 1953!) or mutual information (IScutari and Brogini 
2012|) . Tests in common use do not require any tuning parameter except for the 
type I error threshold a. The estimated B(j) will satisfy © only approximately 
because of type I and type II errors. The former arise from the noisiness inherent 
to the data and limited sample sizes, while the latter are typical of weak depen- 
dencies which will often be omitted from the Markov blanket. 

Several computationally-efficient heuristic algorithms for Markov blanket es - 
timation are available in literatur e, including Grow-Shrink (GS; lMargaritisl 12003b. 



Incre mental Association (IAMBi lTsamardinos et aUl2003l) and Hiton-MB (lAliferis et al. . 
20101) . Conditional independence tests are performed in order of increasing com- 
plexity, thus ensuring that in practice only a small number of SNPs is used for 
each test. 



2.2 Kinship Estimation 

In the past, pedigree information was used to specify kinships, but such informa- 
tion is often missing or inadequate. SNP-based methods for measuring kinships 
have become increasingly common and have the advantage of measuring the re- 
alised amount of genom e sharing, as opposed to t he expected value provided by 
pedigree-based methods (|Astle and Bald ing. l2009h . 



The SNP-based kinship of two individuals is usually based on the average over 
SNPs of the product of their genotypes, coded as 0, 1 and 2 according to the count 
of one of the two alleles. In the following, we denote this genotype matrix with 
X, with rows corresponding to individuals and columns to SNPs, and with X,- its 
z'th column. 
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In human genetic s, kinship is commonly me asured as the proportion of shared 
alleles at each locus (|Morris and CardonL |2007[) . This approach is also known as 
identical-by- state (IBS) kinship, and will be denoted by Kq. Unlike other kinship 
matrices below, Ko is always non-negative and cannot be expressed in the form 
XiXf which leads to parameters directly interpretable as SNP effect sizes (see 
Section [3]for details). 

A nother choice , common in plan t and a nimal genetics, is to centre the geno- 
types ( Habier et al. . 2007 . VanRaden . 2008|) and estimate the kinship matrix as 



i in 

K, =-Y i (Xi-2p i ){Xi-2pi) 



m : , 

i=l 



(2) 



where pi is a vector with every entry equal to the population allele fraction, usually 
estimated as the mean of X{/2. Centring improves interpretability, since kinship 
values can be interpreted as an excess or deficiency of allele sharing compared 
with random allocation of alleles, and so zero can be interpreted as "unrelated". 
However, the requirement to esti mate the pu usually from the same data set, can 
cause problems in some settings (|Astle and BaldingL 12009b . 

One criticism of both the above choices is that the sharing of a rare allele be- 
tween two individuals counts the same as the sharing of a common allele. One 
natural approach to giving more weight to the sharing of a rare allele is to stan- 
dardise over SNPs, thus obtaining 



K 9 



i m 

m ti 



where 



X: 



X,-2p, 



(3) 



The (i, j) entry of K2 can be interpreted as an average over SN Ps of the correlation 
coeffi cient estimated from a single pair of individuals, i and j (lAstle and Balding , 
2009h . 



A modification of K2 has been recently proposed by ISpeed et al.l (|2012al) . 
based on evidence that the effects of SNPs are sensitive to uneven LD across 
the genome. In particular, SNP effects are over-estimated in high-LD regions 
and under-estimated in low-LD regions. To correct for this bias, SNPs can be 
re-weighted: 



(4) 



where the weight vector w 



[wi •• -Wp] solves 



mm 

w 



EH 

!=1 



C,w| 



(5) 
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and C, is a vector of squared correlations of SNP i with neighbouring SNPs. 
Weights decay exponentially with physical distance, according to a decay rate 
X whose value reflects average LD for the data set. For computational reasons, 
the minimisation in © is performed separately on different chromosomes and, 
within each chromosome, on different regions chosen based on A. 



3 Materials and Methods 



We explored the effects of the approaches outlined in Section |2] on the predictive 
power of GS models using three p ublicly-available real- world data sets. The yield 
data from the AGOUEB project dWaugh et all l2010L ICockram et ail l2010h con- 
sist of 227 UK winter bar l ey var i eties and 810 SNPs . The heterogeneous mouse 



population (|Solberg et al.L I2006L IValdar et all |2006|) from the Wellcome Trust 



Case Control Consortium (WTCCC) consists of 1940 SNP profiles and 12545 
SNPs; among the recorded traits , we consider growth rate and weight. The rice 
data set from Izhao et al.l (|2011|) consists of 413 varieties of Oryza sativa with 
73808 SNPs; among the 34 recorded traits, we consider the number of seeds per 
panicle because of its low variability among the various subpopulations included 
in the original analysis. 

All data sets have been preprocessed by removing SNPs with minor allele 
frequencies < 1% and those with > 20% missing data. The missi ng data in the 



remai ning SNPs have been imputed using the impute R package (IHastie et al 



2012J). Other wid ely used imput ation methods in genetics, such as that imple- 
mented in MaCH (Li et all l201Ch . were not available because of the absence of 
accurate SNP maps; the position of many SNPs is unknown, and only genetic dis- 
tances (in cM) were available between mapped SNPs. Furthermore, we removed 
one SNP from each pair with correlation above > 0.9 to increase the numerical 
stability of the models. 

To investigate Markov blanket feature selection, we considered the following 
GS models: 

• Ridge regression, LASSO and the elastic net penalised regressions. These 
are all based on 



y = i u+X/3 + e with j8 



argmin{Ai||/3| 

j8 



i+A 2 ||/3|| 2 },A 1 ,A 2 ^0, (6) 



where y is the trait of interest, X are the SNP genotypes, /3 are the fixed 
SNP effects and £ are independent, normally-distributed errors with vari 



ance Of . We used t he implementation prov ided by the penalized (jGoeman, 



2012|) and glmnet ([Friedman et al.L 1201 Of) R packages. When considering 
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the elastic net we restricted both the L\ and L2 penalties for the genetic ef- 
fects j5 to be strictly positive (Ai, A2 > 0), to facilitate the comparison with 
ridge regression (Ai = 0) and the LASSO (A2 = 0). 

Partial least squares (PLS) regression as implemented in the pis R package 
jMeviket all 120 lib - 



Gene tic BLUP (GBLUP) implemented in the synbreed R package (|Wimmer et al 
20121) . It uses the linear mixed model 



y = i u+Zg + £, 



g~7V(0,Kc7, 



(7) 



where g are the random effects and Z is a design matrix that can be used for 
example to indicate the same genotype exposed to different environments 
(IPiepho et all 120121) . Any positive definite matrix can be used for K. 

When K can be expressed in the form XX r , GBLUP can be shown to be 
equivalent to the Bayesian linear regression 



y = £ x* j8j + £ with SNP effect prior # 




(8) 



in which K determines the transformation X* of the SNP genotypes. For 
instance, the original X are used when K = Ki ; the scaled X, from © when 
K = K2; and the weighted w,X,/£w;- from © when K = K3. This formula- 
tion of GBLUP results in a more natural interpretation of SNP effects, and 
is sometimes known as random regression BLUP (RR-BLUP). 

Markov blanket feature selection has been p erformed with the IAMB algo- 
rithm as implemented in the bnlearn R package (IScutarill2010|) . using the exact 
Student's t test for Pearson's correlation with a type I error threshold of a = 0. 15. 
Each GS model was fitted both using all the available SNPs and using only the 
SNPs included in the Markov blanket. 

The different kinship matrices were investigated within GBLUP, as the other 
GS models do not include an explicit kinship term. Kj and K2 were c omputed us 



ing syn breed. For K3, we used the freely available LDAK software (|Speed et al 
2012bh . The LD decay rate was set to A = 50cM for the AGOUEB data, A = 
0.2cM for the mouse data and A = lOOcM for the rice data. Such values were 
found, through experimentation, to reduce the bias introduced by LD without af- 
fecting the g enetic information p resent in the SNP profiles. Ko was computed 
with PLINK dPurcell et all 120070 . All configurations of GS models and kinships 
were fitted once using all available SNPs and once using only the Markov blanket. 
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The predictive power of the GS models was measured with Pearson's correla- 
tion coefficient p between the observed trait values and the predictions obtained 
from 10-fold cross-validation. For each model, cross-validation was run 5 times 
and the resulting correlations averaged. The correlation between observed and 
fitted trait values is also reported as a measure of goodness of fit. 



4 Results 

Tabled] reports the observed correlations (p, i.e. the correlation between the ob- 
served and the fitted trait values) and the predictive correlations (pcv> i- e - me cor- 
relations obtained from cross-validation) for PLS, ridge regression, LASSO and 
the elastic net. The corresponding correlations arising from the subset of SNPs 
included in the Markov blankets are labelled pMB an d Pcvmb, respectively. 

First of all, we note that for a = 0.15 Markov blankets only select a small 
number of SNPs, regardless of the dimension of the SNP profile. The average size 
of the Markov blankets obtained from cross-validation is 185 for the AGOUEB 
data, 543 (for growth rate) and 525 (for weight) for the mouse data, and 293 for 
the rice data. Of those SNPs, 136 (74%) appear in at least half of the cross- 
validation folds for AGOUEB, 241 (46%) for the mouse data and weight, 276 
(51%) for the mouse data and growth rate, but only 15 (5%) for the rice data. 
This can be attributed to the very low ratio between sample size and number of 
SNPs in the rice data (< 0.01) compared to the mouse (0.15) and AGOUEB (0.28) 
data. As expected, the dimension reduction is smaller in the case of the AGOUEB 
data because of the limited num ber of available SNPs , despite the extensiv e LD 



present in cultivated UK barley (|Cockram et all I2010L iRostoks et all |2006J). On 



the other hand, only a small proportion of the original SNPs are retained for the 
mouse and rice data sets (about 4% and 0.4%, respectively). In each case, the 
number of SNPs included in the Markov blankets is smaller than the sample size, 
thus ensuring the regularity and numerical stability of the GS models. 

The position of mapped SNPs within the respective genomes is shown in Fig- 
ure \T\ (AGOUEB), Figure |2] (rice), Figure |3] (mice, weight) and Figure |4] (mice, 
growth). For all but the AGOUEB data, we can clearly see how the Markov blan- 
kets arising from cross-validation identify some regions as associated with the trait 
being modelled (e.g. SNPs with mapped in the range [0.1cM,75.8cM] of chro- 
mosome 1 are included with high probability for both traits in the mice data set) 
while completely discarding other regions (e.g. [96.3cM, 108. 3cM] in chromo- 
some 2 and [70.6cM, 88.5cM] in chromosome 3). The positions of these regions 
may provide useful prior information in subsequent association studies and in tar- 
geting future geno typing efforts. In the case of the AGOUEB data, marker density 
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Model 


P 


Pmb 


A 


Pcv 


PCVMB 


A 


AGOUEB, YIELD (185 SNPs out of 810, 23%) 


PLS 


0.812 


0.805 


-0.007 


0.495 


0.495 


+0.000 


Ridge 


0.817 


0.765 


-0.051 


0.501 


0.489 


-0.012 


LASSO 


0.829 


0.811 


-0.018 


0.400 


0.399 


-0.001 


Elastic Net 


0.806 


0.752 


-0.054 


0.500 


0.489 


-0.011 


MICE, GROWTH RATE (543 SNPs out of 12545K, 4%) 


PLS 


0.716 


0.882 


+0.166 


0.344 


0.388 


+0.044 


Ridge 


0.841 


0.889 


+0.047 


0.366 


0.394 


+0.028 


LASSO 


0.717 


0.881 


+0.164 


0.390 


0.394 


+0.004 


Elastic Net 


0.751 


0.893 


+0.142 


0.403 


0.401 


-0.001 


MICE, WEIGHT (525 SNPs out of 12545K, 4%) 


PLS 


0.927 


0.823 


-0.104 


0.502 


0.524 


+0.022 


Ridge 


0.877 


0.843 


-0.034 


0.526 


0.542 


+0.016 


LASSO 


0.743 


0.807 


+0.064 


0.579 


0.577 


-0.001 


Elastic Net 


0.789 


0.845 


+0.056 


0.580 


0.580 


+0.000 


RICE, SEEDS PER PANICLE (293 SNPs out of 73808K, 0.4%) 


PLS 


0.853 


0.923 


+0.070 


0.583 


0.601 


+0.018 


Ridge 


0.950 


0.921 


-0.029 


0.601 


0.612 


+0.011 


LASSO 


0.885 


0.939 


+0.054 


0.516 


0.580 


+0.064 


Elastic Net 


0.958 


0.917 


+0.040 


0.602 


0.612 


+0.010 



Table 1: Correlation coefficients for PLS, ridge regression, LASSO and the elas- 
tic net: p is the correlation between observed and fitted trait values; pcv is the 
predictive correlation obtained from cross-validation; p MB and Pcvmb are the cor- 
responding quantities obtained using only the SNPs in the Markov blanket. The 
highest value for each quantity and data set is shown in bold. The average di- 
mension of the Markov blanket over cross-validation is reported in parentheses 
for each data set and trait. 
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Figure 1: Frequency of the SNPs included in the Markov blankets learned from 
the AGOUEB data using cross-validation, plotted against the position of the SNPs 
in the barley genome. Green ticks indicate the positions of all mapped SNPs for 
this data set. 
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Figure 2: Frequency of the SNPs included in the Markov blankets learned from 
the rice data using cross-validation, plotted against the position of the SNPs in the 
genome. Green ticks indicate the positions of all mapped SNPs for this data set. 
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Figure 3: Frequency of the SNPs included in the Markov blankets learned from 
the mouse weight data using cross-validation, plotted against the position of the 
SNPs in the barley genome. Green ticks indicate the positions of all mapped SNPs 
for this data set. 
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Figure 4: Frequency of the SNPs included in the Markov blankets learned from 
the mouse growth data using cross-validation, plotted against the position of the 
SNPs in the barley genome. Green ticks indicate the positions of all mapped SNPs 
for this data set. 
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AGOUEB 


MICE, 


MICE, 


RICE 


Kinship 






GROWTH 


WEIGHT 






matrix 


P 


Pcv 


P 


Pcv 


P 


Pcv 


P 


Pcv 


K 


0.848 


0.511 


0.838 


0.376 


0.931 


0.536 


0.933 


0.596 


Ki 


0.847 


0.512 


0.656 


0.366 


0.882 


0.507 


0.933 


0.590 


K 2 


0.848 


0.513 


0.688 


0.388 


0.883 


0.508 


0.933 


0.598 


K 3 


0.832 


0.521 


0.695 


0.400 


0.881 


0.554 


0.918 


0.594 



Table 2: Correlation coefficients obtained in GBLUP using the four kinship ma- 
trices defined in Section 12.21 The highest value for each quantity and data set is 
shown in bold, p and pcv are defined as in Tabled] 



is not high enough to identify regions with markedly different association levels. 

We observe no loss in the predictive power of the GS models after the Markov 
blanket feature selection. In fact, the increased numerical stability resulting from 
the reduced number of SNPs slightly improved the predictive power of the S mod- 
els. The average of pcv over the four analyses was 0.481, 0.498, 0.471 and 0.521 
for PLS, ridge, LASSO and elastic net respectively, while the corresponding av- 
erages for pcv,MB were 0.502, 0.509, 0.487 and 0.520, all with an approximate 
standard deviation of 0.0057 computed as in lHoopen (|1958l) . 

If we choose a < 0. 15, we obtain Markov blankets that are too small to capture 
polygenic effects (results not shown). A possible explanation for this behaviour 
may be that large values of a allow Markov blankets to initially include SNPs that 
are weakly associated with the trait, to the point that they would be individually 
discarded. In addition, among them there may be sets of SNPs that are jointly 
significant due to epistasis, and such sets are retained in the Markov blanket. 

Furthermore, Markov blankets outperform other subsamples of the same size. 
To show this, we generated for each data set 100 random subsets of SNPs of the 
same size as the corresponding Markov blanket. Subsequently, we used them to 
fit the GS models and to compute the corresponding Pcvmb- As we can see from 
Figure |5l the Markov blanket always results in higher values of pcv,MB man me 
random subsets. 

The elastic net consistently outperforms the other GS models both with and 
without the use of Markov blankets, except for the AGOUEB data set (in which 
pcv is essentially the same for ridge regression and the elastic net). 

Overall, from Table |2] we see that the predictive performance of GBLUP im- 
proves as the kinship matrices progress from Ki through to K3. Ko, while not 
being competitive with K3, outperforms at least one of Ki and K2 for all data 
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AGOUEB MICE, WEIGHT MICE, GROWTH RICE 




0.1 0.2 0.3 0.4 0.5 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.3 0.4 0.5 0.6 




0.1 0.2 0.3 0.4 0.5 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.3 0.4 0.5 0.6 




0.1 0.2 0.3 0.4 0.5 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.3 0.4 0.5 0.6 




0.1 0.2 0.3 0.4 0.5 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.3 0.4 0.5 0.6 



Figure 5: Comparison between the cross-validated correlations obtained from the 
Markov blankets (pcv.MB, vertical red dashed line in each panel) and those ob- 
tained from randomly selected subsets of SNPs of the same size. 
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\{~\ ncnin 


PCV,MB 


A 


P MB, KIN 


A 


PCV,MB 


A 


P MB, KIN 


A 


matn y 


AGOUEB, YIELD 


RICE, SEEDS/PANICLE 


^0 


0.418 


-0.093 


0.479 


-0.032 


0.426 


-0.170 


0.597 


+0.001 




0.412 


-0.100 


0.482 


-0.030 


0.428 


-0.161 


0.592 


+0.002 




0.414 


-0.099 


0.491 


-0.022 


0.429 


-0.168 


0.589 


-0.008 


K 3 


0.415 


-0.105 


0.475 


-0.045 


0.425 


0.169 


0.592 


-0.003 




MICE, GROWTH RATE 


MICE, WEIGHT 


K 


0.194 


-0.182 


0.378 


+0.002 


0.219 


-0.317 


0.534 


-0.002 


Ki 


0.118 


-0.248 


0.357 


-0.008 


0.120 


-0.387 


0.457 


-0.005 


K 2 


0.176 


-0.211 


0.363 


-0.025 


0.182 


-0.326 


0.480 


-0.028 


K 3 


0.195 


-0.204 


0.379 


-0.021 


0.225 


-0.328 


0.530 


-0.024 



Table 3: Correlation coefficients for GBLUP using the four kinship matrices de- 
fined in Section 12.21 and Markov blanket feature selection. Pmb an d Pcvmb are 
defined as in Table [Q Pmb,kin is the predictive correlation obtained from cross- 
validation with the use of Markov blankets but with the kinship matrices estimated 
from the full SNP profile. The highest value for each quantity and data set is 
shown in bold. Differences are computed from p (for Pmb) and pcv (for pcv,MB 
and Pmb, kin) in Table H 



sets but AGOUEB. The means of the five p C v values are 0.504 for K , 0.493 
for Ki, 0.501 for K2, and 0.518 for K3, all with an approximate standard devi- 
ation of 0.0057. Thus, GBLUP with K3 performs as well as the elastic net and 
outperforms PLS, ridge regression and the LASSO. 

Although the elastic net performed equally well with or without Markov blan- 
ket feature selection, that is not the case for GBLUP (Table |3]). For all kinship 
matrices, the reduced size of the Markov blanket relative to the full SNP set de- 
tracts from the computation of kinship coefficients, leading to a substantial loss of 
predictive power. If all SNPs are available and can be used to compute the kinship 
matrices, then much but not all of this loss is restored. 

5 Conclusions 

We have shown that Markov blanket feature selection applied as a preliminary step 
in GS is able to greatly reduce the size of the SNP set with no loss (and possibly 
a small gain) in the predictive power of PLS, ridge regression, LASSO and the 
elastic net. Among those models, the elastic net was the best performer, followed 
by ridge regression. If GS is to be performed repeatedly for the same phenotype, 
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for example in successive generations of crops, Markov blanket feature selection 
opens the possibility of reducing costs by genotyping many fewer markers. 

In the absen c e of a feature selection step, the LD-adjusted kinship matrix 



K3 (|Speed et all l2012af) provide s slightly better predictiv e power than the ma- 



trix with no LD adjustment K2 ( Astle and Balding . 2009 ) and the IBS kinship 
matrix Ko produced by PLINK (IPurcell et all 120070 . In turn, K2 and Ko appear 
superior to the matr ix with neither LD adjustment nor standardising of SNPs K] 
(jHabier et all [2007b- Using K3, GBLUP was competitive with the elastic net 
(both had mean pcv = 0.52 over the four datasets). 

Markov blanket feature selection is not compatible with GBLUP because of 
the requirement for large numbers of SNPs to compute the kinship matrix. How- 
ever, Markov blanket feature selection has only a small adverse effect on GBLUP 
if all SNPs are available for computing the kinship matrix. 
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